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Emerets’s experiments on pressurized sulfur hydride have shown that H 3 S metal has the highest 
known superconducting critical temperature Tc = 203K. The Emerets data show pressure induced 
changes of the isotope coefficient between 0.25 and 0.5, in disagreement with Eliashberg theory 
which predicts a nearly constant isotope coefficient. We assign the pressure dependent isotope 
coefficient to Lifshitz transitions induced by pressure and zero point lattice fluctuations. It is 
known that pressure could induce changes of the topology of the Fermi surface, called Lifshitz 
transitions, but were neglected in previous papers on the H 3 S superconductivity issue. Here we 
propose that H 3 S is a multi-gap superconductor with a hrst condensate in the BCS regime (in the 
large Fermi surface with high Fermi energy) which coexists with a second condensates in the BCS- 
BEC crossover regime (located on a small Fermi surface spots with small Fermi energy) near the F 
and M point. We discuss the need of Bianconi-Perali-Valletta (BPV) superconductivity theory for 
superconductivity in H 3 S. It includes both the correction of the chemical potential due to pairing 
and the configuration interaction between different condensates, neglected by the Eliashberg theory. 

Here the shape resonance in superconducting gaps, similar to Feshbach resonance in ultracold gases, 
gives a relevant contribution to amplify the critical temperature. Therefore this work provides some 
key tools needed in the search for new room temperature superconductors. 

PACS numbers: 74.20.Pq,74.72.-h,74.25.Jb 


I. INTRODUCTION. 

Following early claims of 190 K superconductivity in 
sulfur hydride at very high pressure [l|, new results of 
near room temperature superconductivity with Tc=203 
K i.e., at only -70 °C 0 have been presented on June 17, 
2015 at Superstripes 2015 conference in Ischia, Italy 0 
and they have triggered a very high scientific interest . 
The recent work of Emerets’s group [3] shows the Meiss¬ 
ner effect and the pressure dependent critical tempera¬ 
ture of H 3 S and D 3 S. These results have triggered today 
the materials research for room temperature supercon¬ 
ductors in different hydrides at extreme high pressures 

[ft. 

The experimental discovery of near room temperature 
superconductivity in sulfur hydride H 2 S at very high 
pressure was predicted by to occur in the high pres¬ 
sure metallic H 3 S phase, with Imim lattice symmetry. 
Disproportion from 2 (H 2 S) + H 2 to 2 (H 3 S) occurs at very 
high pressure. The theoretical prediction of Duan et al. 
(l 0 l| has been obtained by using the successful theoretical 
approach used to predict crystalline structures at high 
pressure in material science: the evolutionary algorithm 
Universal Structure Predictor: Evolutionary Xtallogra- 
phy (USPEX) [HI- 

Many of preceding theoretical studies on this issue con¬ 
clude that the superconducting phase in pressurized H 3 S 
is described by the Eliashberg theory [TO, [TMl while 
Hirsh proposes the hole superconductivity model [20l| . 

The Eremets’s group research was motivated by the 
search for room temperature superconductivity predicted 


to emerge in metallic hydrogen and hydrides [2TI - l29l |. 

The BCS theory [sO] has given a microscopic descrip¬ 
tion of the superconducting condensate many body wave- 
function made of interacting Cooper pairs in a weak cou¬ 
pling regime, where the pairing is mediated by the con¬ 
ventional attractive phonon-exchange mechanism. We 
call here standard BCS theory the BCS formulas [s^ 
obtained with many approximations, valid for a simple 
homogeneous crystal, with a single band and isotropic 
pairing, using a single value of the density of states at 
the Fermi level Nq and a constant coupling constant A. 
Moreover the standard BCS theory assumes a very small 
energy of the phonon and a very high electron density 
i.e., a high Fermi energy oj^/Ep << 1 called Migdal ap¬ 
proximation [sij l within the adiabatic Born-Oppenheimer 
approximation, where the electronic and ionic degrees of 
freedom can be rigorously separated. The prediction of 
the superconducting critical temperature in the frame of 
the standard BCS theory approximations has required 
the introduction of the electron-electron repulsive inter¬ 
action and corrections due to strong cou plin g given by 
the McMillan [^ [s^ and the Eliashberg |34| formula. 

It was rapidly well accepted that the critical tempera¬ 
ture cannot be larger than 30 K in the frame of standard 
BCS theory based only on the role of high energy phonons 
and strong electron-phonon coupling [H, [T^, . In fact 

in the single band approximation increases with both 
phonon energy and coupling strength but for extreme 
strong electron-phonon coupling the electron liquid at 
low temperature prefers to order in the real space, form¬ 
ing electronic crystals like, charge density waves, spin 
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FIG. 1: (Color online) The crystalline structure of H 3 S with 
ImSm lattice symmetry, with two formula units per unit cell, 
with sulfur (yellow large spheres) and hydrogen atoms (red 
small spheres). The linear S-H-S hydrogen bonds along the 
a,b and c axis in the /m3m lattice of H 3 S, form a first 3D 
network of three linear chains of covalent bonds crossing at 
the sulfur atom at S1(0,0,0,) which coexists with a second 3D 
network of other three linear chains of covalent bonds (solid 
black lines) crossing at the sulfur atom at 81(0.5,0.5,0.5). 


density waves, which compete with the superconducting 
phase. Moreover if superconductivity survives increasing 
electron-phonon coupling, the critical temperature de¬ 
creases since the phonon energy is pushed toward zero. 
In these regime the lattice structure collapses and the 
system is in the verge of a catastrophe. Therefore it was 
known that materials research for room temperate su¬ 
perconductivity cannot be driven simply by looking to 
increase the phonon energy and coupling strength in the 
cooper pairing. The theoretical predictions of high 
in solid hydrogen and hydrates at high pressure were 
based not only on the high frequency phonons mediat¬ 
ing the pairing in solid hydrogen or hydrates but also on 
joint control of the electron electron interaction via the 
changes of the dielectric constant controlling the repul¬ 
sive electron-electron interaction and the Coulomb en¬ 
ergy [ 2 ll - [^ which can become negative m, [H m for 
an electronic system in the low density limit. 

The self consistent quantum many body theory of su¬ 
perconductivity, avoi ding the BCS approximations, was 
made by Bogoliubov [35l| . giving the fundamental spec¬ 
trum of excited quasiparticles. It was developed by 
Gorkov and Blatt [13] including the contribution to 
the critical temperature of the condensation energy, re¬ 
lated with Josephson-like [i^ terms, [H, and consid¬ 
ering multiple symmetries of possible multiple conden¬ 
sates with a single critical temperature [dlj. A second 


route to rise the critical temperature of superconductiv¬ 
ity toward the maximum possible energy, KbTc/E p ^ 1 
was proposed in 1994 based on the non standard BCS 
theories where the control of can be achieved not only 
by increasing the pairing strength but also by the control 
of the condensation energy, via a pair exchange mecha¬ 
nism between condensates, called the shape resonance 
in superconducting gaps [4a - [45l| . It was proposed that 
cuprates are actually complex inhomogeneous materials 
with lattice quenched disorder, lattice modulations and 
short range charge density waves giving multiple Fermi 
surface arcs where in each Fermi arc there is a different 
superconducting condensate. One of these condensates 
is in the polaronic regime, where the Migdal approxima¬ 
tion is violated ojq/Ep ~ 1 , and the condensate is in the 
BCS-BEC crossover regime while the other condensates 
are within the Migdal approximation loq/Ep << 1. In 
1996-1998 [iol - lsoj a superconductivity theory was devel¬ 
oped by Bianconi, Pera^ Valletta (BPV) based on the 
Blatt and Legget theories. The superconduct¬ 
ing critical temperature here is not only controlled only 
by the cooper pair formation but also by the exchange 
terms between pairs. These terms control the condensa¬ 
tion energy and phase coherence and are unavoidable in¬ 
gredients for the formation and stability of high temper¬ 
ature superconducting condensate with short coherence 
length. A similar mechanism, called Feshbach resonance, 
was proposed later to increase the critical temperature 
for the formation of the superfluid condensate in ultra¬ 
cold atoms Later other many body theories have 

been proposed aimed to control the condensation energy 
and the global phase coherence in superconductors and 
superfluids [S^Im]. 

Details of the normal metallic phase become key ingre¬ 
dients in these non standard BCS theories: a) the com¬ 
plex Fermiology, beyond the single band model; b) the 
formation of charge density waves and polaronic Wigner 
crystals involving the electronic component in n-th Fermi 
surfaces beyond the Migdal approximation uio/Enp ^ 1; 
c) the strong electronic correlations in the electron fluid 
with very low Fermi energy, beyond the usually assumed 
Fermi gas approximation; e) the anomalous or negative 
dielectric response controlling the electron-electron in¬ 
teraction; f) the complex inhomogeneous spatial geom¬ 
etry in systems with nanoscale phase separation |55l | in¬ 
duced W electron-lattice interaction and lattice misfit 
strain |56l | which can giv e also insulator to superconduc¬ 
tor phase transitions [57j |. 


II. SUPERCONDUCTIVITY IN BINARY 
INTERMETALLICS: A15 AND DIBORIDES 

The overall features of superconductivity in H 3 S show 
that it is like other binary intermetallics, as diborides 
and A15 systems therefore we shown in hgure 2 the evo¬ 
lution of the maximum critical temperature of diatomic 
intermetallics. 
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FIG. 2: (Color online) The superconducting critical temper¬ 
ature of superconducting elements, open black circles, binary 
intermetallics, filled blue triangles, and ternary and quater¬ 
nary oxide layered perovskites made of Cu02-atomic layers 
with different X spacer layers (filled red squares) coefficient 
as a function of the year of the discovery. 


A15 compounds (like NbsGe) have the same ImZm lat¬ 
tice structure as H 3 S above 100 GPa (shown in Fig.l) 
which is made of two intertwined networks of atomic 
wires running in 3D as noticed by Friedel (H^ . within 
a 3D density of states (oTI-l^ and superconductivity ap¬ 
pears in highly inhomogeneous phases at the edge of a 
lattice catastrophe due incipient structural phase transi¬ 
tions (63 - 1^ giving coinplex pattern of defects and local 
lattice fluctuations [671 - 1^ . 

Following the discovery of superconductivity in MgB 2 
in 2001 it was first proposed the breakdown of the 
Eliashberg theory and the need of BPV theory including 
shape resonance to describe this system made of mul¬ 
tiple condensates and a single critical temperature. The 
breakdown of the Eliashberg theory was rapidly accepted 
by the community. The theoretical assumptions used by 
standard BGS theory, considering a single band and sin¬ 
gle gap, fail due to anisotropic pairing in the clean limit. 
The standard BGS formulas failed to predict accurately 
the unusually high transition temperature, the effects of 
isotope substitution on the critical transition tempera¬ 
ture, and the anomalous specific heat of MgB 2 . While 
some authors proposed the multi-band superconductiv¬ 
ity theory in the clean limit and weak coupling, it 
was rapidly shown the need of strong coupling anisotropic 
two band Eliashberg theory to describe this unexpected 
novel superconducting phase [zl, [t^ . The isostructural 
AIB 2 system with only tt electrons at the Fermi level 
becomes superconductor only when the chemical doping 
pushes the top of the a band above the chemical poten¬ 


tial giving a Lifshitz transition for the appearing of two 
small hole-like a Fermi pocket s 173 . Moreover because of 
zero point atomic fluctuations [75l [^. the Fermi energy 
in the a bands is time and space dependent therefore the 
Migdal approximation ujq/Ep < 1, is violated and the 
^ Eliashberg theory fails, when the Fermi energy is tuned 
^ over a large e nergy r ange of 600 meV above the bottom 
of band edge.jo, 13 Therefore in spite of the conven¬ 
tional phonon mediated pairing, the 40 K superconduct¬ 
ing phase in magnesium diboride needs a non standard 
BGS theory. Gonsidering all data collected in diborides 
doped with Sc and A1 for Mg or with C for B it was 
possible to give a theoretical interpretation of the vari¬ 
ation of the two superconducting gaps and the critical 
temperature as a function of the energy separation be¬ 
tween the top of the a band and the chemical potential 
using the BPV theory [t^. The superconducting a gaps 
are much larger than the gaps in the large tt Fermi sur¬ 
faces and evolve as a function of the energy separation 
of the top of the a band and the chemical potential!^. 
The small percentage of the partial electronic density of 
states, DOS(cr), relative to the total DOS(fof), usually 
entering in the standard McMillan formula for T^, gives a 
large superconducting a gap, which drives the full system 
to the highest critical temperature known in binary inter¬ 
metallics, before the recent discovery of Eremets’s group 
Q. In a single band metal, where oJojEp ~ 1 the Migdal 
and the adiabatic approximation s 17611 f ail, the system 
enters in the BGS-BEG crossover and both weak 

coupling BGS [13 and strong coupling Eliashberg theo¬ 
ries breakdown. On the contrary the superconductivity 
phase in a system with multiple condensates, where only 
one or two of the Fermi surface spots are in the BGS-BEG 
ctossover while other condensates in the large Fermi sur¬ 
faces are in the BGS regime, can be described by the 
theory of Bianconi Perali Valletta (BPV) [Hj, l79j cor¬ 
rectly describing the multi-gap superconductivity at the 
BGS-BEG crossover avoiding standard approximations. 
The condensation energy of the pairs is determined by 
many body configuration interaction between pairs form¬ 
ing the single macroscopic superconduc ting quantum co¬ 
herent phase. The condensation energy [3 is relevant a) 
in the BGS-BEG crossover An/Epn 1 and in the for¬ 
mation of the superconducting quantum coherent phase 
with a single Tc in a system made of multiple gaps A„ in 
different bands with different Fermi energies Epn- The 
shift of the chemical potential going from the normal to 
the superconducting phase below T^, becomes a key term 
in MBBC, while on the contrary it is considered to be 
negligible in standard BGS. The results of the BPV the¬ 
ory applied to doped magnesium diboride [zll show that 
it is in a regime of multiple condensates with different 
symmetry with the key role of the shape resonance in 
the superconducting gaps between one condensate in the 
BGS-BEG crossover, in the a band, and other conden¬ 
sates in the BGS regime, in the tt bands, with an essential 
Josephson-like pair exchange term. 
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FIG. 3: (Color online) The pressure dependent isotope coef¬ 
ficient as a function of pressure (filled circles) calculated by 
interpolation of the experimental data of the critical temper¬ 
ature of HsS (dashed green line) and D 3 S) solid blue line 
reported by Droz et al. [31 


III. LIFSHITZ TRANSITIONS. 


The topological changes of the Fermi surface of the nor¬ 
mal phase caused by the lattice strain (as a response to 
pressure), misfit strain, chemical substitution, variation 
of the electronic density) are called Lifshitz transitions 
[s^ . In many cases Lifshitz transitions appear on points 
or lines of high symmetry in the Brillouin zone and are as¬ 
sociated with changes of the symmetry or dimensionality 
of the wave functions of the electrons at the Fermi level 
in critical Fermi surfaces. The Lifshitz transitions are re¬ 
vealed experimentally by anomalies in lattice parameters, 
in the density of states near the Fermi energy, in elastic 
properties, in anomalous thermodynamic and transport 
properties of metallic materials [83|. Lifshitz [s^. noted 
that at zero temperature, T = 0, Lifshitz transitions are 
true phase transitions of order 2.5 (as in Ehrenfest’s clas¬ 
sification) therefore are called ”2.5 Lifshitz transition”. 
A sharp Lifshitz transition at T=0, at high temperature 
shows a crossover character. In presence of strong inter¬ 
actions the 2.5 phase transition becomes first order with 
a phase separation between two phases where the chem¬ 
ical potential in each phase is shifted above and below 
the Lifshitz transition respectively[^[83,ll3 driving the 
system well tuned at a Lifshitz transition on the verge of 
a lattice catastrophe as shown in the case of A15 and di¬ 
borides. While in the early times the interest was focused 
on Lifshitz transitions in single band metals, now the in¬ 
terest is addressed to Lifshitz transitions in multi-band 
metals with different Fermi surfaces showing multi-gap 


superconductivity. The shape resonance mechanism con¬ 
siders the relevant contribution, near the Lifshitz transi¬ 
tion, of the exchange interaction between pairs in the hot 
spots and in all other points of the Fermi surfaces in the k 
space. We have learn that the ’’Devil is in the details”, in 
fact for the optimization of a high critical temperature it 
is necessary to reach particular Lifshitz transitions where 
a small number of electrons in a new appearing Fermi 
surface (the hot spot) are in the extreme strong coupling 
regime in the antiadiabatic uiQ/Epn ^ 1 regime in the 
normal phase without the lattice catastrophe since the 
majority of the electron gas is in large Fermi surfaces 
well in the Midgal approximation uj^/Epn < 1. The 
shape resonance in superconducting gaps is a Josephson- 
like term describing a contact interaction between pairs 
which increases the critical temperature for the pairs con¬ 
densation, very similar to the Feshbach resonance in ul¬ 
tracold gases driying up the ratio between k^Tc 

and the Fermi energy up to value of 0.2. In the super¬ 
fluid phase in the hot spot is in the BCS-BEC crossover 
An/Epn ^ 1 [ 1 ^ and the shape resonances in the super¬ 
conducting gaps can be optimized to amplify the critical 
temperature to the highest value (50j | . A large zero point 
lattice motion is a key term in the theory of the shape res¬ 
onance in superconducting gaps whichj is in action in the 
proximity of the Lifshitz transitions. It was proposed to 
play a key role in cuprates 4 ^ ^ ISTl - fs^ . it was veri- 


lied be in action in diborides 7417911 a nd it was confirmed 
in iron based superconductors [90l - l^ . We discuss below 
the failure of the standard Midgal approximation and 
of the breakdown of Eliasberg theory in H 3 S. The large 
electron-phonon coupling and the high frequencies of H- 
phonons contribute to rise Tc [E lisl . IE but we claim 
the need of theories beyond Eliasberg theory, like the 
multigap BCS-BEC crossover (MBBC) theory and the 
BPV shape resonance mechanism to describe the gaps 
in the new small Fermi surfaces appearing at the Lifshitz 
transitions controlled by pressure and their configuration 
interaction with the gap in the large Fermi surface. 


IV. ISOTOPE EFFECT IN H3S. 

The isotope effect in H 3 S [l| has provided a direct ev¬ 
idence of the involvement of the lattice degree of free¬ 
dom in the pairing process. Therefore the isotope effect 
in sulfur hydrate at high pressure H 3 S has been inter¬ 
preted as ruling out theories of unconventional supercon¬ 
ductivity (based only on spin liquid models or magnetic 
interactions) and supporting conventional theories of su¬ 
perconductivity based on the role of lattice fluctuations. 
However the standard BCS theory predicts a pressure in¬ 
dependent isotope coefficient 0.5 while in non standard 
BCS theories, like in the multigap anisotropic BCS, the 
isotope coefficient devia tes fr om 0.5 as in magnesium di¬ 
boride where it is 0.26 tl00|. From the new results re¬ 
ported by Droz et al. [2| we have extracted the isotope 
coefficient as a function of pressure shown in figure |3| 
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We can see in figure [3] large variations of the isotope 
coefficient reaching a minimum of 0.2 and first maxi¬ 
mum reaching 0.5 at 170 GPa and a second peak at 240 
GPa. The anomalous pressure dependent isotope coef¬ 
ficient has been found in cuprates superconductors as a 
function of doping [H, with anomalies where the 
chemical potential crosses Lifshitz transitions driven by 
pressure. Therefore the data in figure |3] indicate the pos¬ 
sible presence of Lifshitz transitions in the pressure range 
showing near room temperature superconductivity. 


V. BAND STRUCTURE CALCULATION OF 
HsS AS A FUNCTION OF PRESSURE 

We hav e pe rformed preliminary band structure 
calculations 1 10 ij of H 3 S with Jm3m lattice symme¬ 
try mad e using the linear muffin-tin orbital (LMTO) 
method fTn^ . lIo3l and the local spin-density approxima¬ 
tion (LSDA) |l04j| . We show in Fig. |T]the dispersion of 
the bands crossing the chemical potential in the T — M 
direction. 

Self-consistent paramagnetic calculations are made for 
a simple cubic unit cell containing 8 sites totally, used 
for A15 compounds with the same crystalline structure. 
The deta ils of the method have been published earlier 
fioMo^ . The present calculations are in goo d agree - 
ment with previous band structure calculations [I5l - ll9l| . 
We confirm the presence of a narrow peak of the occupied 
total Density of States (DOS) very close to the chemical 
potential at 200 GPa i.e., for a=5.6 a.u. The narrow peak 
in the total DOS in a narrow energy range around the 
chemical potential is shown in Fig. [SJ While previous 
papers have stated that this peak is pinned at the zero 
energy, we show in Fig. [5] that the energy position of 
this peak relative to the chemical potential shifts with 
pressure. The narrow peak of the DOS is pushed toward 
high energy by pressure and it crosses the chemical po¬ 
tential at the highest P for the lattice constant smaller 
than a=5.8 a.u.. 

The narrow peak of the DOS peak is related with the 
flat dispersion of bands in the T — M direction in the en¬ 
ergy range of 2 eV below the chemical potential, shown in 
Fig. m The band structure shows a first steep band with 
very large energy dispersion with its band edges at about 
20 eV away from the chemical potential. This first band 
gives the large Fermi surface shown in [l3l| . There are 
three other small Fermi surface pockets centered around 
the F point shown in Fig. |4]for different lattice parameter 
a. 

The tops of the hole-like bands near the F point are be¬ 
low Ep at low pressure P but they move above the chem¬ 
ical potential above 130 GPa pressure as shown in Fig01 
The edges of these bands cross the chemical potential as 
function of P giving topological Lifshitz transitions with 
the appearing of new Fermi surfaces. Moreover there is 
a small Fermi surface hole pocket due to the crossing of 
one band at about 2/3 of the T — M distance indicated 



FIG. 4: (Color online) The dispersion of the bands in the 
electronic structure of H 3 S calculated using the simple cubic 
cell, used for A15 compounds with ImSm structure, between 
M and F points. The band dispersion in this direction gives 
the narrow DOS peak near the Fermi energy. Moreover we 
show the variation of this bands at different pressure with 
changing the lattice parameter between a=6.2 a.u and a=5.6. 
We show that there are bands (d and 7 ) forming 3 small 
Fermi surfaces at the F point and another band crosses the 
zero energy at the point indicated by S which gives a new 
small hole pocket at high pressure. 



FIG. 5: (Color online) The total DOS for H 3 S at different lat¬ 
tice constants. A sharp peak of the DOS crosses the chemical 
potential at high pressure. 
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by S in Fig. 0] which appears only for lattice parameters 
smaller than a ~ 5.8. a.u., i.e. in the pressure range of 
highest P when is highest. The top of this band goes 
from -0.2 eV below the chemical potential to -1-0.1 eV 
when a decreases from 6.2 to 5.6 a.u.. Finally we asso¬ 
ciate the crossing of the chemical potential by this band 
at S in Fig. |4]with the Fermi level crossing of the narrow 
DOS peak in Fig. [S] 


Let us now consider the fact that because of the zero- 
point motion (ZPM) the different Fermi energies in the 
small pockets show large energy fluctuations like in mag¬ 
nesium diboride. The ZPM is large because of the small 
mass of H atoms and it has large effects on the fluc¬ 
tuations of electronic structure. Such effects have been 
shown to be important in several different m aterials, even 
if their atomic masses are larger [IIMIl- Lattice fluc¬ 
tuations can perturb spin waves and phonons in high-T^ 
cuprates and it cannot be neglected if superconductiv- 
ity r elies on few phonons coupled with particular bands 
[ll5l| . The Debye temperature for H phonons is high 
(~ ISOOiti) and the amplitude of lattice fluctuations from 
ZPM is large already at low T. With a force constant 
K = Muj"^ of 7 eV/A^ we obtain an average amplitude u 
of the order 0.15 A. As seen in Fig. [SI the first high ly¬ 
ing valence band has a width is about 2 Ry. This makes 
the band dispersion and Fermi velocities high in the first 
band forming the Fermi surface centered at the R point. 
The effects of energy band broadening is negligible here 
since the chemical potential is far from band edges. The 
large effect of the zero point motion is on the small Fermi 
pockets near the F point. The low-T energy band fluctu¬ 
ations in materials with narrower band widths has been 
fo und to be ab out 20 meV for u in the range 0.03-0.04 
A fiilln^ [n5|. From an extrapolation of these values to 
the conditions in H 3 S we estimate that the band energy 
fluctuation can be of the order of 160 meV for H-bands. 
Therefore when the chemical potential is tuned by pres¬ 
sure near a Lifshitz transition, so that the Enp in one 
of the bands is of the order of 160 meV the topology of 
the small Fermi surfaces made of small hole- or electron- 
pockets shows strong dynamical fluctuations controlled 
by the zero point lattice fluctuations. 


Figure [5] shows the energy position of the top of the 3 
bands at the F point and the band at about 2/3 of the F — 
M distance as a function of pressure. The energy spread 
due to zero point lattice fluctuations is indicated by the 
red area for states at the F point by the yellow area for the 
states at F — M point. The shift of the energy position of 
the narrow peak in the DOS near the chemical potential 
tracks the top of the band at the T-M point indicates 
that the chemical potential here is tuned at the Lifshitz 
transition appearing around 170 GPa for the appearing of 
this new Fermi surface. Finally these results show that in 
superconducting pressurized sulfur hydrate metal is made 
of one large Fermi surface and four small pockets tuned 
by pressure at Lifshitz transitions like doping tune the 
metal in cuprates at Lifshitz transitions for the appearing 
new Fermi surface arcs [1^ . 
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FIG. 6: (Color online) Energy shift of the tops of the three 
hole-like bands at the F point (solid red line, dashed red line 
and dotted yellow line) crossing the chemical potential at 100 
GPa, where the errors bars and the red region indicates the 
zero point energy fluctuations of the band edges. The filled 
black dots and the dotted black line show the energy shift of 
the top of the band crossing the chemical potential at about 
2/3 of the V — M distance of the Brillouin zone. The open 
circles and the blue dashed line show the energy shift of the 
narrow peak of DOS near the chemical potential 


VI. BPV THEORY FOR MULTI-GAPS 
SUPERCONDUCTORS AT THE BCS-BEC 
CROSSOVER 


The standard BCS theory in undergraduate courses 
considers a homogeneous single band metal and describes 
the cooper pairing in weak coupling limit as due to quan¬ 
tum exchange of a low energy phonon between two high 
energy electrons at the Fermi level in the Migdal approx¬ 
imation where the gap energy is much smaller than the 
Fermi energy and the gap to critical temperature ratio is 
equal to 1.75. Non standard BCS theories diverge sub¬ 
stantially from the standard BCS theory. Let us simplify 
a long story by classify 7 different non standard BCS 
theories. Let us start with the first group considering 
a single effective band, in a dirty limit, giving a single 
condensate: lA) the standard strong coupling Eliash- 
berg theories considering a high energy phonon or vi- 
bron (a quantum of intramolecular vibration) interacting 
with electrons having a large Fermi energy loq/Ef << 1. 
Most of previous theories for high Tc in sulfur hydrates 
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have proposed this approach for a large phonon frequency 
and strong coupling; IB) the theories considering the 
case of a low density electron gas beyond the Midgal ap¬ 
proximation loq/E pn > 1 where the single condensate is 
formed by BEC condensation or by BCS-BEC crossover 
An/Ep ~ 1; 1C) non standard BCS theories which con¬ 
sider pairing mediated by electronic excitations, called 
exciton theories. 

The second group of non standard BCS theories consid¬ 
ers the anisotropic BCS theory, in the clean limit, there¬ 
fore they focus on multi-gap superconductivity where 
multiple gaps are formed in different in different spots 
of the k-space. These non standard BCS theories con¬ 
sider metals with multiple bands with different symme¬ 
try crossing the chemical potential: 2A) The multiband 
BCS model where all multiple Fermi surfaces are in the 
BCS limit; 2B) the Fermi-Bose model where the Fermi 
level of a first band is degenerate with a single level, occu¬ 
pied by paired electrons. Here a BCS condensate coexist 
with bosons which undergoes a Bose condensation; 2C) 
the extreme case of superconductivity driven mostly by 
exchange like interband pairing with very weak, or no 
intraband pairing; 2D) The theory of shaperesonance in 
superconducting gaps for multi-condensates with a small 
interband pairing, emerging in metals where: a first small 
charge density in small Fermi surfaces form condensates 
in the BEC- BCS crossover, beyond Migdal approxima¬ 
tion, which coexists with a majority of charges in large 
Fermi surfaces forming BCS condensates. 

A third group of theories consider the case of sys¬ 
tems with relevant electronic and lattice inhomogeneity 
with insulator to superconductor transitions in presence 
of nanoscale phase separation where the Josephson-like 
interaction between localized condensates play a key role. 

Here we focus on the BPV (Bianconi, Perali, Valletta) 
theory which correctly inculde the shape resonance (2D) 
mechanisin, p ro p osed for cuprates, [^. HH liM ISTHSa 
diborides and iron based superconductors [90l - l9^ . 

This mechanism (2D) is proposed here for sulfur hydrides 
at very high pressure following the reported evidence for 
Lifshitz transitions driven by high pressure. 

We recall bellow the BPV theory for multi gap super¬ 
conductors at the BCS-BEC crossover including shape 
resonances. Let us consider a system made of multiple 
bands with index n, having a steep free electron like dis¬ 
persion in the x direction and a flat band-like dispersion 
in the y direction. The energy separation between the 
chemical potential and the bottom of the n-th band de¬ 
fines the Fermi energy of the n-th band. This formula¬ 
tion was proposed for complex systems in presence of 
one-dimensional lattice modulation or one-dimensional 
charge density waves where the chemical potential is 
tuned near a Lifshitz transition like in magnesium di¬ 
borides, A15, cuprates, iron based superconductors and 
we propose here for sulfur hydrate at high pressure. 

The formula for the superconducting critical temper¬ 
ature Tc in a anisotronic multi-gap non standard BCS 
scheme is given by the linearized BCS equation consider¬ 


ing the simplest case of a two dimensional system but 
the extension to three dimensional system is trivial. 
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where the energy dispersion is measured with respect to 
the chemical potential. 

We consider a superconductor with multiple gaps 
^n,ky in multiple bands n with flat band-like dispersion 
in the y direction and steep free-electron-like dispersion 
in the x direction for a simple model of a two dimensional 
metal with a one-dimensional superlattice modulation in 
the y-direction. The self consistent equation for the gaps 
at (T = 0) where each gap depends on the other gaps is 
given by 
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where N is the total number of wave-vectors in the 
discrete summation, /i is the chemical potential, V^'^, is 
the effective pairing interaction 
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Here we take account of the interference effects be¬ 
tween the wave functions of the pairing electrons in the 
different bands, where n and n'are the band indexes, 
ky{k'y) is the superlattice wave-vector and kx{k'^) is the 
component of the wave-vector in the free-electron-like di¬ 
rection of the initial (final) state in the pairing process. 
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Here Nq is the DOS at Ep without the lattice mod¬ 
ulation, An,n' is the dimensionless coupling parameter, 
S = LxLy is the surface of the plane and ipn,ky{y) are 
the eigenfunctions in the ID superlattice. The gap equa¬ 
tion need to be solved iteratively. The anisotropic gaps 
dependent on the band index and on the superlattice 
wave-vector ky. According with Leggett [4l| the ground- 
state BCS wave function corresponds to an ensemble of 
overlapping Cooper pairs at weak coupling (BCS regime) 
and evolves to molecular (non-overlapping) pairs with 
bosonic character and this approach remains valid also if 
a particular band is in the BCS-BEC crossover and be¬ 
yond Migdal approximation because all other bands are 
in the BCS regime and in the Migdal approximation. 







However in this anomalous regime, where Eliashberg 
theory breakdown, on density: by increasing coupling or 
decreasing the density by approaching the band edge, 
the chemical potential fi results strongly renormalized 
with respect to the Fermi energy Ep of the non interact¬ 
ing system, and approaches minus half of the molecular 
binding energy of the corresponding two-body problem 
in the vacuum. In the case of a Lifshitz transition, de¬ 
scribed in this paper, all electrons in the new appear¬ 
ing Fermi surface condense forming a condensate in the 
BCS-BFC crossover. Therefore at any chosen value of 
the charge density for a number of the occupied bands 
Nf,, the chemical potential in the superconducting phase 
should be renormalized by the gap opening at any chosen 
value of the charge density p using the following formula: 
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taking the increment in ky as Sky = 27:jLy for a size 
of the considered surface &Lx Ly and in the range 


f-min = max [o, p-UJo- En^ky] , 
e-max = max [O, p + uio - E„^ky] , 



VII. CONCLUSION. 

In this work we have presented the breakdown of the 
Fliashberg theory for H 3 S, in fact the electronic struc¬ 
ture of sulfur hydrides H 3 S with Imim lattice structure 
as function of pressure shows Lifshitz transitions revealed 
by band crossings at Ep and by the shift of the narrow 
peak in the density of states below the chemical potential 
pushed above it by lattice fluctuations associated with 
the hydrogen zero point motion. We have discussed the 
presence of two topological Lifshitz transitions at two 
critical pressures Pci = 110 Gpa, and Pc 2 = 175 GPa 
by pressure dependent electronic structure calculations 
of H 3 S. At the first Lifshitz transition, around Pci = 110 
Gpa, three new Fermi surface spots appear at the 7 point 
pushed up by pressure. These anisotropic bands are char¬ 
acterized by a flat dispersion in the T — M direction and 
a steep dispersion in the T — R direction. The second 
Lifshitz transition at Pc 2 = 180 GPa is due a new where 
the experimental critical temperature is nearly constant. 


The amplitude of the energy fluctuations of this band 
edge due to atomic zero point motion has been calculated 
and we have found that it pushes this DOS peak above 
the chemical potential. Therefore dynamical energy fluc¬ 
tuations of the band edge due to zero point motion of the 
hydrogen atoms is of high relevance. We find a colossal 
zero point energy fluctuation which induces a 160 meV 
energy fluctuation of the Lifshitz transitions. The present 
results show that the condensates in the 4 small hole 
pockets around the T point and the small Fermi surface 
in the T — M direction are beyond the Migdal approxima¬ 
tion loq/Ep 1 , including lattice dynamical zero point 
fluctuations. In particular the condensate in the hole-like 
6 Fermi surface pocket, associated with the sharp quasi- 
ID peak in the DOS, pushed at the chemical potential 
by pressure, is clearly in the BCS-BFC crossover regime, 
coexisting with other condensates in the BCS regime in 
other large Fermi surfaces in the pressure range where 
the superconducting critical temperature is near room 
temperature. 

The emerging scenario is pairing in sulfur hydrides at 
high pressure in a dynamical landscape where key en¬ 
ergy parameters have all the same magnitude i.e, of the 
order of 160 meV: 1. the energy separation between the 
average position of the n‘^ band edge and the chemical 
potential, dehned as the n*^ Fermi energy controlling the 
’’appearing or disappearing Fermi surface spot” Lifshitz 
transition in the Fermi surface topology, 2. the energy 
separation of a peak in the Density of States and the 
band edge, usually controlling the ’’Neck opening” Lif¬ 
shitz transition 3. the amplitude of the energy fluctu¬ 
ation of the critical Fermi surfaces associated with zero 
point lattice fluctuations 4. The energy of the pairing 
interaction dehning the energy of the cut-off for the for¬ 
mations of pairs away from the chemical potential 

In this scenario the BCS approximations used in the 
standard BCS theory are no more valid and the crit¬ 
ical temperature is controlled not only by the energy 
of the changed boson, wq, and the effective electron- 
phonon coupling (given by the product of the density of 
states times the electron-phonon coupling constant) but 
also by the condensation energy. While in the standard 
Fliashberg theory, the correction to the chemical poten¬ 
tial induced by electron-phonon coupling is ignored on 
the verge of the Lifshitz transition this correction, which 
has much impact, is considered in the Multi-gaps BCS- 
BFC Crossover theory. In this situation the chemical 
shift from the normal to the condensed phase below Tc 
is no more negligible, and the coupling should be renor¬ 
malized by a factor, given by the quantum overlap of 
the condensed pairs in cuprates [4^ |50, M diborides 

[zlliilzi and iron based superconductirs [90l - [^ . 

Finally this paper shows the breakdown of Fliash¬ 
berg theory for pressurized hydrides, supports the role 
of phonons [l^, [l5l - [l^ but the presence of Lifshitz tran¬ 
sitions tuned by pressure need the use of the Multi-gaps 
BCS-BFC Crossover theory including shape resonances. 

Further work is needed to investigate i) the divergent 
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amplitude of lattice fluctuations near the R3m to ImSm 
2nd order structural transition around 180 GPa, ii) the 
large mass difference between H an S which requires the 
consideration of different amplitudes of u for the two 
types of atoms, and it should allow for E and k depen¬ 
dences for the energy fluctuations of Lifshitz transitions. 
The electronic band calculations should be extended to 
large supercells needed for more precise estimates of en¬ 
ergy fluctuations of the electronic structure associated 
with spacial structural fluctuations. We think that the 


discovery of superconductivity in sulfur hydrates very 
near room temperature has narrowed the number of pos¬ 
sible road maps toward new functional superconducting 
materials. Further fundamental research on the mech¬ 
anism of room temperature superconductivity in these 
new phase of matter are needed to clarify this emerging 
physical scenario and they will allow the definition of a 
protocol for the material design of new functional room 
temperature superconductors. 
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